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Abstract 

We study the determination of neutrino mass hierarchy, through neu- 
trino experiments within the next ten years. The T2K neutrino oscilla- 
tion experiment will start in 2009. In the experiment the high intensity 

beam from JHF is directed to Super-Kamiokande (SK) detector 295 
km away. The NOuA (off axies neutrino oscillation) experiment is being 
planned, with the beam from Fermi-Lab directed to a site 6I0km away, 
which is 0,7,14 milliradian off-axies. Both the above experiments will 
measure — > v e oscillation probability. The double-CHOOZ experiment 
under construction detects TT^s emitted by nuclear reactors both through a 
near detector (150m) and a far detector (1.05 Km) and measure X^T — s> zTJ 
survival probability. In this paper, we outline a procedure to determine 
the sign of A31 from the simulated data of the above experiments. 



1 Introduction 

Recent advance in neutrino physics observation mainly of astrophysical obser- 
vation suggested the existence of tiny neutrino mass. The experiment and ob- 
servation have shown evidences for neutrino oscillation. The Solar neutrino 
deficit has long been observed [1,2,3,4), the atmospheric neutrino anomaly has 
been indicates that neutrino are massive and there is mixing in lepton sector 
and currently almost confirmed by KamLAND [8] and hence indicates that neu- 
trino are massive and there is mixing in lepton sector. Since there is mixing 
in lepton sector, this indicate to imagine that there occurs CP violation effect 
in lepton sector. Several physicist have considered whether we can see CP vi- 
olation effects in lepton sector through long baseline oscillation experiments. 
The neutrino oscillation probabilities, in general depends on six parameter two 
independent mass squared difference A21 and A31 , there mixing 612, 623 and 
013 and one CP violating phase 5. There are two large mixing angle (#12, #23) 
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and one small (#i 3 )angle, and two mass square difference, A, 
rriij is the neutrino masses, 
where 



ij = - , with 



A21 = A so i, 



ar, 



(1) 



A31 = AAtmo. 



(2) 



The sign of A 3 iand of 623 when 823 7^0, can not be determine with the 
existing data. For the mass square difference, there are two possibility, A 31 > 
or A31 < 0, correspond to two different types of neutrino mass order normal 
mass hierarchy, mi < m 2 < m 3 (A31 > 0), and inverted hierarchy, mi > m 2 > 
m 3 (A 3 i < 0). The angles 0i2and 823 represent the neutrino mixing angles 
corresponding to solar and atmospheric neutrino oscillation, much progress has 
been made to words, determining the values of the three mixing angle. From 
measurement of the neutrino survival probability i/„ — > v e and v e — > v e in the 
atmospheric flux, so that one mixing angle is near ^and one is small [11] from 
the v e — > v e survival probability in the solar flux, so that the mixing angle is 
either large (LMA) or small (SMA) by solar solution [12]. Nothing is known 
about CP violating phase. In this paper, we tested the sign of A 3 iby using 
three different baseline (T2K, NOvA and Double CHOOZ) experiments. The 
purpose of this paper is to determine the sign of A31 by \ 2 analysis. Section 2 
describes the mixing angle and masses difference. Section 3 describes the mass 
hierarchy effect in — > v e oscillation probability. In Sec-4. Determine the sign 
of A31 by x 2 analysis. Section 5 summarizes the results and conclusions. 

2 Mixing Angles and Neutrino Mass Squared Dif- 
ferences 

The first evidence is the observation of zenith-angle dependence of atmospheric 
neutrino defect [13] dependent of the atmospheric neutrino — > transition 
with the mass difference and the mixing as 



The second evidence is the solar neutrino deficit [14] .Which is consistent 
with — > v T /v e transition. The SNO experiments [15] are consistent with the 
standard solar model [16] and strong suggest the LMA solution. 



Solar neutrino experiments (Super-K, GALLEX, SAGE, SNO and GNO) 
show that neutrino oscillations, neutrino oscillation provide the most elegant 
explanation of all the data [17].. 



A31 = (1 - 2) x l0- 3 eV 2 ,sin 2 28 2 3 = 1.0. 



(3) 



A 2 i = 7 x lO~ 5 eV 2 , sin 2 26 12 = 0.8. 



(4) 



Asoh 



ar 



= 7±5. 3 x 10" 5 ey 2 



(5) 
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tan 2 6 solar = 0.4±°-] 4 . (6) 

Atmospheric neutrino experiments ( Kamiokande, Super-K ) also show that 
neutrino oscillation. The most excellent fit to the all data [17]. 

A atmo = 2.0±i;° 2 x l(T 3 e V 2 , (7) 

sin 2 29 atmo = 0.4t8;& (8) 

The CHOOZ reactor experiment [18] gives the upper bound of the third 
mixing angle #i3as 

sm 2 6»i 3 <0.20 for \ A 31 \ = 2.0 x 10- 3 eV 2 , (9) 

sin 2 6 13 < 0.16 for \A 31 \ = 2.5 x lO^e^ 2 , , (10) 

sm 2 6»i 3 <0.14 for \ A 31 \ = 3.0 x 10- 3 eV 2 , (11) 

at the 90 % CL. The CP phase S has not been constrained. In future neutrino 
experiments, which plan to measure the oscillation parameter precisely. 



3 Mass Hierarchy Effect in Neutrino Oscillation 
Probability 

Let us briefly recall our present knowledge of neutrino oscillation parameters. 
There are three flavors of neutrinos and they mix to form three mass eigensates. 
This mixing is given by 




where mixing matrix U parametrized [19] as 



U = R(6 23 )nR(6 13 )n*R(6 12 ). (13) 

In the above mixing matrix, II is a diagonal matrix containing the CP violat- 
ing phase S and R(6ij)is is the form of rotation matrices. The mass eigenstates 
vi have eigenvalues . Neutrino oscillation probabilities depend on the two 
mass squared differences A 2 i = m| — m 2 ,A 3 i = m\ — m\ , the three mixing 
angles #12, #23, #13 and the CP violating phase S. Solar neutrino data and Kam- 
LAND experiment determine A21 and 612. Atmospheric neutrino data and K2K 
and MINOS experiments determine |A 3 i| and 23 - CHOOZ experiment and so- 
lar neutrino data constrain 0i3to be small. There is no information at present 
on the CP phase 5. The future experiments are expected to measured #13 and 
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determine the sign of A3iand the magnitude of CP violation, in addition to 
improving the precision of known neutrino oscillation parameters. 

Double CHOOZ experiment is a reactor based experiment dedicated to mea- 
suring 13. In this experiment systematic errors are minimised by having iden- 
tical near and far detectors at distances 150 meter and 1050 meter from the 
sources respectively. This experiment can measure non zero value of 6*13 if 
sm 2 2#i3 > 0.05[20]. Daya Bay reactor experiment will have a similar sensi- 
tivity [21]. In T2K experiment the high intensity of beam from JPARC 
accelerator is directed to SK detector 295 Km away. The detector is 2° off-axis 
from the beam, which lead to the neutrino flux peaking at lower energy. T2K is 
a very high statistics experiment that is expected to start taking data in 2009. 
The neutrino flux is about 100 times the flux of K2K. The number of charged 
current events expected, in the case of no oscillation, is about 3100 per year. 
This experiment will improve the precession of A31 and #23 by measuring the 
muon neutrino survival probability P{v^ — > v e ). It can also measure #13 through 
the measurement of P Me - NOvA is also an accelerator based experiment, which 
uses v^ beam from Fermilab 810 Km away. The detecting material in this exper- 
iment is a scintillator which gives it an excellent electron detection capability. 
Thus NOvAcwa. make a precise determination of P^ e ,NOvA which is expected 
to start taking data in 2011, also will be placed at an off-axis location. Because 
of the longer distance the flux NOvA is peaked at higher energy compared to 
that of T2K. Matter term, which is proportional to neutrino energy, causes a 
25% change in P pe whereas the change in P pe of T2K is only about 10% [22]. 
If A31 positive Pancreases, whereas if A31 is negative it decreases. Below we 
describe a procedure by which sign of A31 can be determined using the data 
from Double CHOOZ, T2K and NOvA . We will compute the smallest value of 
#13 for which the sign of A31 can be determined independent of the CP phase 
S. 



4 Mass Hierarchy Effect in P™_Oscillation Prob- 
ability 

• P™ with A 2i = 

Neutrino oscillation probability, — > v e in long base line experiments is modi- 
fied by the propagation of neutrino through the matter of earth's crust [23]. It 
increases the oscillation probability for neutrinos if A 31 is positive and decreases 
it A3iis negative. The reverse is true for anti- neutrinos. Here we consider a 
method of determining the sign of A 3 iusing v^be&ms only. 

In three flavor mixing, v^ — > v e oscillation probability is given by 

P,e = S m 2 23S m 2 2tf 13 W ( L27 ^ 3lL ) , (14) 

where A31 in eV 2 , the baseline L is in Km and the neutrino energy E is 
in GeV. In the above equation, we made the approximation of setting A 2 i = 
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0, which made it independent of 6± 3 and the CP phase S. In long baseline 
experiments, the neutrinos propagate through earth's crust which has constant 
density of about 3gm/cc. The oscillation probability modified by matter effect 
is given by 

(1 97A m T \ 
' E 31 J , (15) 

where 

sin20? 3 = ^*», (16) 

^31 

A™ - ^(A 31 cos20 13 - A) 2 + (A 31 sin26 13 ) 2 (17) 
Here A is the matter term and is given by 

A = 2V2G F NeE = 0.76 x l()- 4 p(in gm/cc)E v (in GeV). (18) 

From the expression of P Me in three flavor oscillations, we can compute the 
magnitude of the terms dependent on A 2 i. It turns out that as the CP violating 
phase 5 varies from — tt to 7r, P^ e changes by about 25 %. Therefore setting 
A21 = is not a good approximation for analyzing matter effects in long baseline 
experiments. 

• P™ with A 2 i + 

Exact expression for P^ e with matter effects is very complicated. The expres- 
sion derived using a perturbation expansion with $13 and a — as a small 
parameters works very well for baselines up to 1000 Km [24, 25]. Carrying out 
the perturbation expansion to second order in the small parameters, the follow- 
ing analytic formula for — > v e is obtained with the assumption of constant 
matter density 

P£ = sin^^^sin^A, - 1)A) 

± a^c ^ 1 3^n2 gl 3^ 12 ^n2^3 sm(A)sm(AiA)gm((1 _ 
A\\l - A 1 ) 

acos6cose 13 sin2e 13 sin29 12 sin29 23 . ... . . . . . 

H -r-7- -7- sin(A)cos(yl 1 A)szn((l - A^A) 

Ai{L — Ai) 

a 2 cos 2 9 2 3sin 2 2d 12 sin 2 (A 1 A) 
A% ' 

where a = A 2 i/A 3 i,A = A 31 L/4E 1 A 1 = 2V2G F NeE/A 31 , G F is the 
Fermic coupling constant and n e is the electron density in earth's crust. We see 
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Figure 1: P^ oscillations probabilities vs E for A 2 i = 2.5 x 10 3 eV 2 ,9\3 — 
8°,L=295 km and L=810 km. The middle line is P™(5 = 0°),the upper line is 
P™(8 = +90°) and lower line is P™(<5 = -90°). 
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Figure 2: P Me oscillations probabilities vs E for -2.5 x lO~ 3 eV 2 ,9 13 = 12°,L=295 
km and L=810 km. The middle line is PJ£(6 = 0°),the upper line is PJ£(5 = 
+90°) and lower line is P™(5 = -90°). 
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Figure 3: P^ e oscillations probabilities vs E for A 2 i = 2.5 x 10~ 3 eV 2 ,8 13 = 
10°,L=295 km and L=810 km. The dashed line is P™(5 = 30°) with A 3i 
positive and dot-dashed line is P™(<5 = 75°) with A 3 i negative. 



the above expression depends on three unknown quantities #13, sign of A 3 iand 
the CP phase 8. 

There are two kinds of degeneracies inherent in the three flavor expression 
for P™. The first one occurs due to following reason. Since #13 is unknown, 
the P™ for positive A 31 and smaller 6* 13 can be essentially the same as the P™ 
for negative A31 and larger #13 This is illustrated in fig(l) and fig(2). Precise 
determination of 6*13 by Double CHOOZ can eliminate this degeneracy. There 
is a further degeneracy involving the CP phase S . At present there is no ex- 
perimental information on this phase. Note that Double CHOOZ is completely 
insensitive to 5. For a given long baseline experiment, it is possible to find two 
values of the CP phase, 5+ and 5~ , such that P™(+A 3 i,c5+) = P™(-A 31 ,6-), 
with all other oscillation parameters, including #13 fixed. This is illustrated in 
fig(3) and fig(4). However, the above degeneracy can occur for only one base- 
line length at a time. In fig(3) P™ for T2K is essentially the same for both 
signs of A 3 ibut P™ can distinguish between the two signs of A31. In fig(4) the 
situation between T2K and is reversed. If we have data from two long baseline 
experiments with different baseline then we can resolve the above degeneracy 
independent of the 5 and sign of A31. 

The following method may be used to test the sign of A31. From the experi- 
ment we will get three piece of data of three different neutrino oscillation exper- 
iment (T2K, NOvA, Double CHOOZ). Ne(T2K), Ne(NOuA) and Ne(Double 
CHOOZ) is the data of three different experiments. We tested, whether the 
hypothesis of positive or negative A31 fit the data better. These number will be 
the function of #13 and S which as yet unknown. We compute P™ numerically 
by diagonalizing the matter dependent mass squared matrix foe each energy 
bin. In next section, we discuss the testing of A 31 sign by \ 2 analysis. 
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Figure 4: P^ e oscillations probabilities vs E for A 2 i = 2.5 x 10~ 3 ey 2 ,#i 3 = 
10°,L=295 km and L=810 km. The dashed line is P™(<5 = -90°) with A 3i 
positive and dot-dashed line is P™(<5 = 90°) with A 3 i negative. 

5 Determine the A31 Sign by x 2 Analysis 

We take the combined data from Double CHOOZ [26], T2K [27] and NOvA 
[28] can resolve the sign of A 3 i. This resolution depends crucially on matter 
effects which in turn depend on #i 3 . If #i 3 is unmeasurable small, it is extremely 
difficult to determine the sign of A 31 . Here we address the question: what is 
the smallest value of #i 3 for which the sign of A 3i can be resolved by the data 
of the above three experiments, independent of the value of S. 

Since there is no data yet from any of these three experiments, we simulate 
data for each experiment. In our calculation we fix the values of the following 
neutrino parameters: A 2i = 8.0 x 10~ 5 eV 2 , 8 12 = 34° and d 2 z = 45°. First 
we take |A 3 i| = 2.5 x lO~ 3 eV 2 . The presently allowed range for 6>i 3 is 0° to 
15° and that for the CP phase 5 is -180°to 180°. We pick the true value for 
#i 3 from its allowed range and similarly for 5.. We call these values 9\ r 3 ue and 
S true . W e take A 3 ito be positive and compute the expected number of events 
in each bin of each experiment for 6^ 3 " e and S true . We smear the computed 
event distributions in energy using the energy resolution functions estimated by 
the respective collaborations. The data obtained after the energy smearing, we 
call to be our simulated data, which consists of 92 data point Np lmu , p = 1, 92. 
.Now we take A 3 i to be negative but keep A 3 i the same. We choose test values 
for 9i3 and 8 which we call #*| s * and <5 test .. With these as inputs, we compute 
theoretical values for the number of events in each bin of each experiment. Thus 
we get 92 theoretical expectations, Np est ,p = 1, 92. We compute % 2 between the 
simulated data and the theoretical values 

92 / N s im _ N th\2 

x \5^\ Off) = ]T ^ (20) 

P =i a P 

In the above discussion, p = 1 , 28 are data of Double CHOOZ, p = 29 to 46 
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are data of T2K and p = 47 to 92 are data of NOvA. a p is the error in A^ im .It 
is the square root of the sum of squares of statistical and systematic errors. In 
calculating the statistical error the background contribution to it is taken into 
account. Following the above procedure we compute a set of ~x 2 (5 test , 9[ 3 st ) for 
all allowed values of 9{ 3 st and 5 test . Since Np m and N* est are calculated using 
different signs of A31, we expect x 2 m Esq. (20) to be large. But, x 2 is a 
function of 9\ 3 st and S test . Because of the parameter degeneracies discussed in 
sec(4) it is possible to have small x 2 for 9\% st ^ 9[ 3 ue and 5 test ^ 5 true even if A 3i 
values have opposite signs in the calculation of N^ 1 " 1 and N^ est . In particular, 
we require 9\ r 3 ue to be large enough such that Double CHOOZ will be able to 
measure its value. If the minimum of X 1 {8 test , ^lf 9 *) is greater than 4.0, then the 
two signs of A31 are distinguishable at 95% CL for the given values of 9{ 3 ue and 
S true . If the minimum x 2 (<5* es S ^{% st ) is l ess than 4.0, then the two signs of A31 
can not be distinguished at 95% CL for the given values of 8\ r 3 ue and S true . We 
repeat the calculation for other values of 9[ r 3 ue and S true . We look for values of 
9[ r 3 ue for which the minimum of x 2 ((5' es * , 9[ 3 st ) is greater than 4.0 for all allowed 
values of S true . The minimum of 9\ r 3 ue for which the above condition is satisfied, 
is the smallest value of #13 for which sign of A31 and hence the neutrino mass 
hierarchy, can be determined irrespective of the value of the CP phase 5. 

We take the Double CHOOZ data [26] is divided into 28 bins. The measure- 
ments of the near detector give us the unoscillated neutrino event rate in each 
bin. The expected measurement in the far detector, for each bin, is given by 



dNf ar dN near (L 



x 



dE v dE„ \ Lf a 



2 



P{(^^)},|A 31 |,0r 3 " e ,L^)}. (21) 



For Double CHOOZ, the expected error in energy measurement is much smaller 
than the bin size. Therefore the energy resolution can be taken to be a Dirac 
delta function. Thus the simulated number of events per bin is given by the 
above equation. 

We see that T2K data [27] is divided into 18 bins. The expected electron 
neutrino event rate, in each bin, is given by 

dN dN 

dt = dt P J^(+^i,6\7 e ^ true ,L T2K ). (22) 

T2K collaboration estimates the error in reconstructing the neutrino energy 
to be 100 MeV. We take the energy resolution function R{E v ,E mea ) to be a 
Gaussian with a = lOOMeV.We obtain the smeared event rate per bin by 

dN e , , 

-tt; 1 sim = / J -j=rR{E v ,E mea )dE v (23) 

^^mea Ojllii/ 

Finally we take the NOvA [28] data is divided into 46 bins, for each of the 
off-axis locations Omrdjmrd and 14mrd. We consider one off-axis location at a 
time. As in the case of T2K, the expected electron neutrino event rate, in each 
bin, is given by 
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Figure 5: Plots of Xmin m 9[™ e — S true plane, Omr off-axis location with 
low energy (left) and medium energy (Right) options for NOvA are assumed. 
IA31I = 2.5 x 10~ 3 eV 2 .The symbol are explained in the text. 



dt = ji/^{+^iA7 e ^ trU ^ L NOvA ). (24) 

Again as in the case of T2K, we obtain the smeared event distribution by 
means of a Gaussian resolution function with a = lOOMeV. We assume that 
both T2K and NOvA will run only in neutrino mode for five years. In com- 
puting the numbers for NOvA, we consider various different possibilities: Low 
energy beam with various different off-axis angles and also medium energy beam 
with various different off-axis angles. The theoretical expectation values are 
calculated using eq. (21), eq. (22) and eq. (24) with -A 3 i,^| s * and 5 test as 
neutrino parameters. Note that no smearing is done in calculating theoretical 
expectation values for event numbers. 



6 Summary 

In this paper, we have studied the neutrino mass hierarchy. Our results are 
displayed in fig. (5), fig. (6) and fig. (7). In each figure we give a plot 
of X^ in in f?*™ 6 — 5 true plane. The star symbol represents Xmin < 4-0 square 
represents 4.0 < Xmin < 9.0 the triangle represents 9.0 < Xmin < 16.0 and circle 
represent Xmin > 16-0. In each figure the left panel is generated assuming that 
NOvAwiW run in the low energy option and right panel is generated assuming 
high energy option. Fig. (5) corresponds to Omrd off-axis location of NOvA, 
fig. (7) corresponds to Imrd off-axis location and fig. (7) corresponds to lAmrd 
off-axis location. From the x 2 analysis, we calculate the minimum value of #13 
for which the sign of A31 can be resolved at 95 %CL. For |A 3 i| = 2.5 x 10~ z eV 2 
the low energy option with Omrd and Imrd off axis location seem to have the 
best resolving ability. We repeated our calculation for other allowed value of 
|A3i|.In table 4, we compute the minimum value of 6\ r 3 ue for which the sign of 
A31 could be resolved at 95% CL, independent of the CP phase. We consider 
the Omrd, Imrd and lAmrd off-axis angles of NOvA for different values of A31 
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Figure 6: Plots of Xmin m ^i™ 6 — <5*™ e plane, 7mr off-axis location with 
low energy (left) and medium energy (Right) options for NOvA are assumed. 
|A 3 i| = 2.5 x 10 _3 ey 2 .The symbol are explained in the text. 




Figure 7: Plots of Xmin m ^i™ 6 — S true plane, 14mr off-axis location with 
low energy (left) and medium energy (Right) options for NOvA are assumed. 
IA31I = 2.5 x !CP 3 eV 2 .The symbol are explained in the text. 



IA31I 


Minimum value of (9*™ e 


1.5 x lir^eV' 1 


15° 


2.0 x lO^eV 2 


9° 


2.5 x \0-' A eV 2 


7° 


3.0 x \Q- :i eV 2 


4° 


3.5 x lO~' A eV 2 


4° 


4.0 x IQ-^eV 2 


4° 



Table 1: Minimum value of 6\ r 3 ue , for which the sign of A31 could be resolved 
at 95% CL, independent of CP phase. 
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From the table we see that the minimum value of #i 3 ue for which the sign 
of A 3 iCOuld be resolved at 95% CL independent of the CP phase. This min- 
imum is the same for Omrd and Imrd off-axis angles of the low energy 
option of NOvA. The results are a little worse for the medium energy option of 
NOvA. Determining the type of neutrino mass hierarchy, whether normal or in- 
verted, constitutes one of the fundamental question in neutrino physics. Future 
long baseline experiments aim at addressing this fundamental issue, but suffer 
typically from degeneracies with other neutrino parameters, namely #i 3 and 5 
The presence of such degeneracies limit the sensitivity to the type of hierarchy. 
Many earlier studies focused on the determination of the sign of A 31 by using 
the data of neutrinos and anti-neutrinos from more then one experiment [29, 
30, 31, 32, 33]. In this present paper, we study the possibility of solving the 
neutrino mass hierarchy using only neutrino data of long baseline experiments 
T2K and NOvA and data from Double CHOOZ. We determined, for each al- 
lowed value of |A 31 |, the minimum value of # 13 for which the sign of A 31 could 
be resolved independent of the value of the CP phase. If |A 3 i| = 0.0025eT^ 2 , 
we can rule out the wrong neutrino mass hierarchy at 95 % CL, for the whole 
range 5 true = -180° - 180°, if 9{ r 3 ue > 7.0°. For larger values of |A 3 i| is less 
then 0.002eT^ 2 the neutrino mass hierarchy can not be resolved by the data of 
the above three experiments for any of the allowed values of 6* 13 . 
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